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Caecilians of the family Caeciliidae Rafinesque, 1814 
comprise two genera, Caecilia Linnaeus, 1758, with 33 
species, and Oscaecilia Taylor, 1968, with nine species 
(Wilkinson et al. 2011). They occur in the eastern parts 
of the Andes, from Costa Rica in Central America to the 
south of Bolivia in South America (Maciel & Hoogmoed 
2011, Frost 2015). Both genera occur in Brazil, represent¬ 
ed by Caecilia armata Dunn, 1942, C. gracilis Shaw, 1802, 
C. marcusi Wake, 1985, C. tentaculata Linnaeus, 1758, and 
Oscaecilia hypereumeces Taylor, 1968 (Maciel & Hoog¬ 
moed 2011, SBH 2014). Caecilia, with the largest number 
of species, is distributed from Panama to southern Boliv¬ 
ia in areas with rainforest (Frost 2015). The type species 
is Caecilia tentaculata Linnaeus, 1758, described already 
in the ‘Systema Naturae’ with ‘America as its type locality. 
According to Nussbaum & Hoogmoed (1979), the holo- 
type was lost and two new specimens from Surinam (with¬ 
out specific locality data) were referred to as syntypes by 
Dunn (1942). However, according to Duellman (2012), 
apparently the type was recently located. This species is 
found in terrestrial habitats and is categorized as ‘Least 
Concern’ by the International Union for Conservation of 
Nature (IUCN 2016). Until now, Caecilia tentaculata has 
been regarded as having a distribution linked to the Ama¬ 
zon forest habitat, but here we expand the known distri¬ 


bution to two isolated ‘brejos-de-altitude’, relict mountain 
forests surrounded by dry Caatinga scrubland in north¬ 
eastern Brazil. Isolated ‘brejos-de-altitude’ have some¬ 
times been considered a biogeographical link between the 
coastal Atlantic forest and Amazonian forests. Finding 
C. tentaculata in the area might therefore suggest that it 
is likely present also in Atlantic forest habitats. Herein we 
confirm the species’ identification using both morphologi¬ 
cal and molecular characters and use environmental mod¬ 
elling to identify if the environmental conditions of this 
region are suitable to promote the distribution of C. ten¬ 
taculata. 

In the state of Ceara, five separate ‘brejos-de-altitude’ oc¬ 
cur, including the Massif of Baturite (4°05’-4°4o’ S, 38°3o’- 
39°io’ W) and the Maranguape Mountain (3°54’-4°03’ S, 
38°32’-38°4o’ W), where our specimens were found. They 
are separated by a gap of about 45 km. Four specimens 
were collected at both these sites (Fig. 1, Appendix 1) dur¬ 
ing herpetological surveys, anaesthetised and subsequent¬ 
ly killed by injections with lidocaine hydrochloride, fixed 
with 90% ethanol or 10% formol, and deposited in the 
Collection of Herpetology of the Universidade Federal of 
Ceara (CHUFC). Specimens were identified according to 
Taylor (1968) and Maciel & Hoogmoed (2011) by assess¬ 
ing diagnostic morphological characters. 
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To further confirm their specific identity, and to resolve 
genetic diversity between samples from widely separated 
geographic regions, two specimens, CHUFC A6775 from 
Maranguape Mountain and CHUFC A6543 from the Mas¬ 
sif of Baturite, were sequenced for a region of the mito¬ 
chondrial 12S rRNA gene by means of Polymerase Chain 
Reaction (PCR) amplification, using the primers nSa 
(CTGGGATTAGATACCCCACTAT) / 12SI1 (GAGGGT- 
GACGGGGCGGTGTGT) from Kocher et al. (1989). 
DNA from tissue samples (muscle or liver) was extract¬ 
ed using standard high-salt protocols (Sambrook et al. 
1989). The PCR protocol was based on a 20-pi reaction mix 
containing 4.0 pi of Blend MasterMix (Hot Fire), 0.5 pi 
of 10 pmol primers each, 2.0 pi of DNA, and distilled wa¬ 
ter making up the balance, following Harris et al. (1998). 
The reaction mix was amplified using the following cycles: 
94°C (13 min), 35 x [94°C (30 sec), 47°C (45 sec) and 72°C 
(1 min)], and 72°C (10 min). Positive and negative checks 
were run with each reaction. Both positive products were 
sent for sequencing to Beckman Coulter Genomics. The 
new sequences were submitted to the GenBank database 
under accession number KY745902. 

The potential model of the C. tentaculata distribution 
was assessed using sampling points in databases (Global 
Biodiversity Information Facility - Gbif 2016, Iucn 2016) 
and literature (Dunn 1942, Maciel & Hoogmoed 2011), 
using Maxent software for species habitat modeling. Fif¬ 
ty-two points of non-redundant occurrence were plotted 
across South America, with a precision of io’ (18.6 x 18.6 
= 344 km 2 ), and posteriorly five non-correlated environ¬ 
mental variables were projected on the mesh and used as 
predictors of the distribution of the studied species based 
on its ecology (mean annual temperature, temperature sea¬ 
sonality, precipitation, precipitation seasonality, and pre¬ 
cipitation over the warmest quarter). The environmental 
variables used were downloaded from WorldClim (2016). 


We also constructed a correlation matrix between all varia¬ 
bles and selected only those that were not highly correlated 
(r < 0.9). The performance of the niche model was eval¬ 
uated according to the curve between specificity (lack of 
commission errors) and sensitivity (lack of omission error) 
and AUC, that is, the point that minimizes the relationship 
between sensitivity and specificity of the model, assessing 
their ability to correctly predict the occurrence of the spe¬ 
cies. The AUC value ranges from o to 1, with values > 0.75 
being considered significant (Elith et al. 2006, Phillips 
& Dudik 2008). This analysis was performed with R plat¬ 
form (R Core Team 2016). 

Borges-Nojosa (2006) recorded two specimens of 
Gymnophiona from the Massif of Baturite, but due to 
their poor state of preservation, listed these erroneously 
as Chthonerpeton aff. arii, a more aquatic caecilian species 
usually known from open areas (Cascon & Lima-Verde 
1994). Diagnostic morphological characters (following 
Taylor 1968, Maciel & Hoogmoed 2011) of the newly 
collected specimens clearly indicated that these specimens 
represented Caecilia tentaculata, a terrestrial forest spe¬ 
cies: in particular the visible eyes, the uniform grey or blue 
natural colouration, total body length (up to 850.0 mm), 
eye-tentacle distance (ET) greater than the distance be¬ 
tween nostril-tentacle (NT), total length 30 times larger 
than the body diameter, and numbers of primary annuli 
(PA) (112-131) and secondary grooves (SG) (12-54). A 368- 
bp fragment of the 12S rRNA gene was successfully se¬ 
quenced from each of the two specimens (CHUFC A6543 
and A6775). These proved identical for this fragment, but 
showed five differences (all transition changes) from the 
sequence of C. tentaculata from GenBank (KF540146) 
from Kaw in French Guiana (San Mauro et al. 2014). Such 
a low level of differentiation (1.4%) is typical of intraspe¬ 
cific variation; the next most similar sequence (C. gracilis, 
KF540147) differed by 9%. Limited genetic diversity across 



Figure 1. Specimens of Caecilia tentaculata: CHUFC A6773, live specimen showing the natural coloration (A); and CHUFC A6774, 
lateral (B), dorsal (C), and ventral (D) views of the head, and cloacal opening (E). Photos: D. M. Borges-Nojosa. 
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Table 1. Morphometric (in mm) and meristic data of Caecilia tentaculata from Ceara state. Legend: TL - total length; BW - body width 
at mid-body; HL - head length measured between tip of snout and posterior edge of skull; HW - head width; ET - distance between 
eye and tentacular aperture; NT - distance between nostril and tentacle; PMT - premaxillary-maxillary teeth; PPT - prevomerine- 
palatine teeth; DT - dentary teeth; ST - splenial (inner mandibular) teeth; PA - primary annuli; SG - secondary grooves (= primary 
annuli bearing secondary grooves); Notes: CHUFC A6477 - immature female; CHUFC A6543 - incomplete specimen. 


Character 

Massif of Baturite 

A6477 A6478 

Locality / Voucher Number 

Maranguape Mountain 
A6479 A2601 A6773 

A6774 

A6775 

Eyes 

visible 

visible 

visible 

visible 

visible 

visible 

visible 

Coloration 

uniform 

uniform 

uniform 

uniform 

uniform 

uniform 

uniform 

TL 

286.0 

268.1 

200.3 

405.0 

363.4 

362.6 

321.8 

BW 

9.1 

5.6 

6.5 

- 

9.6 

9.6 

9.8 

Proportion TL/BW 

31.4 

47.7 

30.7 

- 

37.9 

37.9 

33.0 

HL 

9.5 

10.3 

8.5 

11.9 

12.7 

11.6 

11.1 

HW 

7.0 

6.6 

6 

7.8 

8.57 

8.3 

9.3 

ET 

3.5 

4.1 

4.0 

4.1 

4.3 

4.5 

5.0 

NT 

1.8 

1.2 

1.8 

1.3 

1.9 

2.3 

2.0 

PMT 

11 

21 

19 

15 

15 

16 

20 

PPT 

14 

15 

13 

17 

13 

16 

20 

DT 

15 

19 

17 

16 

16 

17 

15 

ST 

2 

3 

4 

3 

4 

3 

4 

PA 

125 

134 

121 

131 

123 

122 

114 

SG 

26 

23 

15 

22 

20 

16 

20 


large geographic distances has been reported in another 
but not more closely related caecilian, Ichthyophis bombay- 
ensis Taylor, i960, from India (Gower et al. 2007) and 
was suggested to indicate recent genetic continuity possi¬ 
bly associated with variables such as a relatively large clutch 
size. However, as ichthyophiids have an aquatic larval stage 
(San Mauro et al. 2014, Kupfer et al. 2016) that will like¬ 
ly contribute to the dispersal whereas caeciliids are direct- 
developers lacking a larval stage, dispersal processes in 
C. tentaculata are different. Detailed comparisons between 
C. tentaculata from across its entire distribution range may 
shed further light on how environmental or other variables 
might be linked to genetic continuity. Maciel & Hoog- 
moed (2011) defined three OTUs (Operational Taxonom¬ 
ic Units) with head length and distance between nostrils 
as the most discriminant characters in C. tentaculata, al¬ 
though these were based on limited differentiation and 
had only 60% correctly classified. Based on the few sam¬ 
ples available from the “brejos-de-altitude” it is not possi¬ 
ble to infer the OTU. 

The new record extend the known species range by 
about 850 km to the east of previously known localities 
(Viseu, State of Para; Nova Vida, State of Maranhao) (Ma¬ 
ciel & Hoogmoed 2011) (Fig. 2), and add another typical 
Amazon forest species to the herpetofauna of the brejos- 
de-altitude in Ceara, similar to Amphisbaena anomala 
(Barbour, 1914) (Borges-Nojosa & Caramaschi 2003). 
These records, in association with other species that have 
a distribution corridor Amazon-'brejos’-Atlantic, such as 
the snakes Lachesis muta (Linnaeus, 1766) (Borges-No¬ 


josa & Lima-Verde 1999), Sibon nebulatus (Linnaeus, 
1758), Imantodes cenchoa (Linnaeus, 1758) (Nascimento 
& Lima-Verde 1989), Drymoluber dichrous (Peters, 1863) 
(Borges-Nojosa & Lima 2001), and the lizard Cercosaura 
0. ocellata Wagler, 1830 (Borges-Nojosa & Caramaschi 
2003), support the hypothesis of the Amazon and Atlan¬ 
tic Forest being biogeographically linked. Interestingly, al¬ 
though there is a very small region corresponding to the 
Brejos, which the model identifies as being appropriate 
(p > 0.6), the environmental conditions required by the 
species currently seem to be missing from the Atlantic For¬ 
est (Fig. 2). However, this does not by default invalidate the 
hypothesis that the species was widely distributed in past 
periods such as the Pleistocene during which climatic con¬ 
ditions facilitated rainforests to expand. The model for the 
current distribution of the species had a good efficiency 
rating (AUC = 0.823). The variables that contributed most 
to it are ‘seasonality temperature’ (40.6%) and ‘precipita¬ 
tion in the warmest quarter’ (47%). 

The ‘brejos-de-altitude’ are humid, elevated regions that 
are isolated within the dry region of Caatinga scrubland. 
With altitudes ranging between 400 and 1,000 m, they are 
typically covered with relict forest vegetation from adjacent 
humid forests (Andrade-Lima 1982, Prance 1982, San¬ 
tos et al. 2006) with different climatic conditions. They 
contain a unique mix of flora and fauna of three different 
biogeographic regions, Amazonian, Atlantic forest, and dry 
Caatinga (Andrade-Lima 1966, Vanzolini 1981, Santos 
et al 2006, Borges-Nojosa & Caramaschi 2003). Various 
endemic amphibians and reptiles have been reported from 
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Figure 2. Map of ecological niche modelling for Caecilia tentaculata. Squares, triangles and stars are Operational Taxonomic Units 
(Maciel & Hoogmoed 2011), circles other localities (GBIF 2016), and numbers are new records (1 - Maranguape Mountain; 2 - 
Massif of Baturite; specimens of CHUFC) used for modelling. The Maxent default output depicts probability values ranging from 0 
to 1, where higher values indicate higher environmental suitability and therefore a higher probability of species occurrence (Phillips 
& Dudik 2008). 


here, almost all of which are categorized as ‘Threatened’ 
as per IUCN criteria due to the relatively small expanse of 
the ‘brejos’ and ongoing environmental destruction. These 
include amphibians such as Adelophryne baturitensis and 
A. maranguapensis Hoogmoed, Borges & Cascon, 1994 
(VU), Rhinella casconi Roberto, Brito & Thome, 2014, 
the lizards Copeoglossum arajara (Rebou^as-Spieker, 
1981), Leposoma baturitensis Rodrigues & Borges, 1997 
(EN), and the snake Atractus ronnie Passos, Fernandes 
& Borges-Noiosa, 2007 (EN) (MMA 2014). The areas of 
‘brejos’ in the state of Ceara are considered ‘areas of highest 
priority’, as they represent an Atlantic Forest relict (MMA 
2000, Brasil 2004). Our new records confirm that the 
‘brejos’ contain a mix of florae and faunae of both Ama¬ 
zonian and Atlantic origins. Finding C. tentaculata, about 
850 km from the nearest known population, highlights 
not only the uniqueness of the ‘brejos’, but also that new 
and detailed herpetological surveys are urgently needed to 
comprehensively understand the ancient link between the 
Amazon and the Atlantic forests, two of the most biodi- 
verse regions of the world. 
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Appendix 1 

Caecilia tentaculata (8 specimens) Brazil: Ceara: Serra de Ma- 
ranguape, Maranguape Municipality: CHUFC A2601 (s/nc); 
CHUFC A6773 (bo86 ); CHUFC A6774 (15284); CHUFC A6775 
(B285), Maci^o de Baturite, Pacoti Municipality: CHUFC 
A6477 (DMB916); CHUFC A6478 (DMB483); CHUFC A6479 
(DMB684); CHUFC A6543 (BAT03). 
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